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In this paper the entropy of the two-channel Kondo impurity is calculated. We demonstrate that 
the exactly solvable system of finite size can be described in terms of free fermions. None of these 
fermions is Majorana fermion. This allows us to find the impurity entropy S(T, L) for arbitrary 
values of T and L. As found before, S(T = 0, L = oo) equals to (1/2) In 2. For finite L the entropy 
shows a cross-over to the infinite-size behavior at T > Tl = Vf/L, in agreement with results of 
other researchers. 
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I. INTRODUCTION 

Recent investigations of macroscopic quantum phe- 
nomenon, ID quantum wires and Kondo effect discover 
that the zero-temperature entropy of some systems pos- 
sesses an anomaly |T]~Q . Specifically, zero-temperature 
entropy of the two-channel Kondo impurity was found to 
be: 



S(T = 0) = iln2. 



(1) 



This result is in obvious contradiction with the fun- 
damental fact: for any system which exists in nature it 
is assumed that exp(5(T = 0)) is the degeneracy of the 
ground state, i.e. integer number. Equation ([[]), how- 
ever, implies that the degeneracy equals to \/2. 

It was pointed out later in the context of Kondo model 
[^],[| that the entropy value (1/2) In 2 is obtained in the 
thermodynamics limit where system size L is sent to in- 
finity. For a finite size system a non-zero temperature 
scale can be defined in the following way: Tl = vp/L. 
Thermodynamic limit calculations are not valid below 
Tl- Mathematically, this situation can be described as 
follows ||: 



arbitrary values of T and L. As a function of T and L the 
entropy demonstrates a cross-over at T ~ Tl, in agree- 
ment with the theory of critical phenomena. For T > Tl 
finite size effects are unimportant while below Tl the im- 
purity decouples from the Fermi sea. We confirm this by 
both numerical and analytical calculations. 

The paper is organized as follows. In Section II we 
calculate the spectrum of the finite size version of Kondo 
model. As we said this model can be reduced to a set of 
free fermions. In Section III the thermodynamics of the 
model is discussed. Section IV contains the conclusions. 



II. SPECTRUM OF TWO-CHANNEL KONDO 
IMPURITY 

In this part we are going to diagonalizc the hamiltonian 
of two-channel Kondo impurity placed in the box of size 
L. 

In paper I two-channel Kondo hamiltonian 



2 pL 

H = iv F ^2 / dxvjj a (x)d x ip ia (x) + Hi 



(3) 
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Only the first of these two limits is physically relevant |7]] 
and it gives reasonable result. 

As shown in || (further cited as I) two-channel Kondo 
model is exactly solvable at particular values of coupling 
constants. In this paper we study that exactly solvable 
model in the limit of finite L. A description of the fi- 
nite size exactly solvable system in terms of free Dirac 
fermions is found. Majorana fermion which appears in 
L = oo system combines with another Majorana fermion 
and forms Dirac fermion. In L = oo limit that second 
Majorana fermion remains hidden in the continuous spec- 
trum. Numerical calculation of the spectrum of the free 
fermions allows us to obtain the impurity entropy S for 
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Jata^vL (0)^(0) (4) 



was studied in the limit of infinite L. Here tp - left- 
moving electrons, a, (3 are electron spin indices, index i 
is electron "flavor", t x is the impurity pseudospin. Hj 
describes the coupling of the electron spin to the impurity 
pseudospin. Coupling constants J x .y,z are chosen in such 
a way that J x — J v , J z = ttvf > 0. 

Following the procedure described in I (bosonization 
and re-fermionization) this hamiltonian can be trans- 
formed to: 



H„ 



iv F / dx^(x)d x ^(x) + 



(5) 



1 



[**(()) + tf(0)] {<# -d), 



where ^Sf is some new fermionic field, d is a fermionic 
operator representing the impurity, a is a smallest length 
in the system (unit cell size). Expressing this hamiltonian 
in terms of c k s, Fourier components of '5, we obtain: 



L new — / 
k 



A + -d), ( 6 ) 



A = 



J. x 



(7) 



Hamiltonian (^) written in terms of Fourier components 
(^) will be a starting point for our calculations. 
To proceed further let us introduce two operators: 



D 1 = -L(it + d), 
D, = ^-d). 



(8) 
(9) 



Obviously, D\ and D 2 satisfy the following relations: 

D{ 2 = Di, a , (10) 
{D 1 ;D 1 } = {D 2 ;D 2 } = 1, (11) 
{£>i;D 2 } = 0, (12) 

Though ( [L2|) i s usual property of fermionic operators the 
equalities (|l0|) and ( |l l| ) define Majorana fermion opera- 
tors. It was pointed out in I that D\ commutes with the 
hamiltonian: 



[H new ;Di] — 0. 



(13) 



The anomalous entropy behavior in case of infinite L was 
attributed to the presence of this Majorana fermion com- 
muting with (JsJ) . We will see below that there is no Ma- 
jorana fermion and no anomalous entropy for the finite 
size system. 

Now let us diagonalize the hamiltonian (^J) . Since it is 
quadratic in fermionic operators, we will look for a set of 
eigenoperators of the form: 



IkCk 



5D 2 , 



such that 



(14) 



(15) 



Conditions ( fjr| - |l5|) give rise to an eigenvalue problem 
for 7's and 8: 



u% = e-klk - V2A5, 
W7fe = -tklk - V2AS, 
jS = -V2A^(7fc+ 7fc ). 



(16) 
(17) 
(18) 



Here e k = vpk. 

The solution of this set of equations is: 



7fe 
Ik 



= 2A 2 V — 



.V2A 

u-e k 

. \/2A 
)— . 

lu + e k 

1 



uj + e k uj - e k 



(19) 
(20) 
(21) 



The condition {/•; /} — 1 fixes the value of S. 

The solution shown above allows one to verify the 
following property of symmetry: equations (|l6|-[l"8|) are 
invariant under the transformation uj — > —uj,^ k ~ > 
~lk,lk —> ~7fc,<5 — * S. Another way of expressing this 
symmetry is: 

fl = f-u>- (22) 

This implies that uj > solutions of equations (nq- 
|l8| ) constitute fermionic spectrum of the hamiltonian (]q) 
while uj < solutions are mirror image of that spectrum. 
In other words, u> > eigenvectors give us fermionic cre- 
ation operators /t while uj < eigenvectors give us de- 
struction operators: f_\ w i = f\u\- 

There is another remarkable property of equations ( [l6| - 
|l^): uj = solution always exists. This can be readily 
seen from (|2l|) . The following is true for ui — eigenop- 
erator (see (1221) and (151)): 



(23) 
(24) 



fu=Q /oj=Q) 

[H new ; fu, = o] = 0. 



Equations (p3|) , ( 53) show us that fu,=o is the Majorana 
operator (see properties ( |To| ) and (11)) commuting with 
the hamiltonian (ph. Having two Majorana operators we 
can define two other operators: 



F = -^(i/«=o + I>i) ) 



(25) 
(26) 



These two operators are annihilation and creation oper- 
ators for some new fermion: 



= 0,c^0. 



(27) 
(28) 



Since both f u —o and £>i commute with hamiltonian, the 
new fermion has zero energy: 



[Hnew] F] — 0. 



(29) 



Thus, the set of eigenoperators of H new consists of f u for 
uj > 0, F and their hermitian conjugates. All of them 
are usual fermionic operators. We have a description of 
our system in terms of free fermions with no Majorana 
fermion. The presence of the zero-energy fermion F im- 
plies that S(T = 0, L) = In 2. 
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III. FREE ENERGY OF THE KONDO IMPURITY 

In order to calculate the entropy let us consider free 
energy of our system. In the previous section we have 
shown that our system can be described in terms of free 
fermions. Therefore, the impurity free energy is: 



T - 



-Tin 2 



T£ln 

iu>0 



1 



j/T 



■To. 



(30) 



where Tq is the free energy without the impurity. The 
impurity entropy is: 



b ~ Of' 



(31) 



We solve (^l|) numerically and find the entropy for dif- 
ferent T. Three curves for three different values of L are 
plotted on fig.l: L a < Lf, < L c = oo. One can see from 
the graph that the system has a cross-over at Tj,: as T 
goes below Tl the impurity decouples from the Fermi sea 
due to the presence of the finite size gaps in the excitation 
spectrum of the conducting electrons. It is clearly seen 
on fig.l that the limiting value of S(T, L) at T = and 
L = oo depends on the path in (T, L) plane, in agreement 
with (0). 




FIG. 1. The impurity entropy vs. temperature (in the 
units of r) for different values of L: L a < Lb < L c — oo, 
L a /Lb = 0.6. 

There is one peculiarity on the graphs of the entropy 
for the finite systems. The entropy of the impurity de- 
creases as a function of temperature below Tl- There 
is no contradiction here with the fundamental criterion 
of thermodynamical stability: the impurity entropy is 
not the entropy of any physical system but rather a dif- 
ference between the entropy of the system with Kondo 
impurity and the entropy of free electrons without impu- 
rity. For both of these two systems the entropy increases 
with temperature, but their difference does not. 



We see form fig.l that the result for L = oo works fine 
as high temperature approximation. In order to show this 
analytically let us calculate high temperature expression 
for the impurity free energy (|30|). 

In the finite size limit the spectrum consists of a set 
of eigenvalues defined by (pij) . The characteristic dis- 
tance between two neighboring eigenvalues is Tl = Vp/L. 
However, for T 3> Tl one can expect that discreteness 
of the spectrum does not play any role and it can be 
approximated by continuum. 

In order to study continuous version of our model we 
are going to obtain the density of states p(uu). To do so 
we will consider the following auxiliary function 



h{uj) 



to 



lu -2A 2 



1 



1 



k 0J + ek 



LO — Ck 



(32) 



JJ(cj + e k )(uu - e k ). 

k 

As it follows from ( pl| ) this function has two important 
properties: if uj belongs to the spectrum of our model 
then h{uj) — 0; h{uS) does not have any singularities. It is 
tempting to wright that p(u>) — (l/7r)(d/dw)Im(ln/i(w)). 
However, this is incorrect. There is one thing we have to 
take care of: for a given eigenvalue lu ^ the function h 
has two zeros (at u> and — u>), for to = it has a second 
order zero. The following is the correct expression for the 
density of states: 



p = — — Im [ — \nh(io 
7T dui \ 2 



(33) 



The factor of (1/2) removes the impact of extra zeros of 
h(u). 
Then: 



1 d 

p = — Im (lnw- 

1 2-Kduj v 



(34) 



hi 



UJ 



2A 



2 E- 
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uj + t k W - £fe 



+ e*:) + J2 ln ( w - efe ) ) • 

k k J 

Expression in square brackets is evaluated in the follow- 
ing way: 



2A 2 V — 



1 



A 2 



2L 



lu + e k uj — e k 



de de 



(35) 



tvvf J-n uj + e uj — e 



lu — iT, 
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where bandwidth f2 and coupling T are defined as follows: 

n = ^, (36) 

a 



P 

r = — 2- < a 



(37) 



After straightforward calculations the density of states 
is found to be 



p{u) = p + 



1 



1 



2tt to 2 + r 2 



8(u), 



(38) 



where p is the density of states without the impurity. 
Thus, the impurity free energy is: 



T 



-T 



-n 



- — 5 — : =5- In ( 1 + e~ u/T \ dw-— In 2. 
27rw 2 + r 2 V J 2 



(39) 



This free energy coincides with the free energy found in 
I up to constant. From ( |3g| ) zero temperature value of 
the impurity entropy is found to be (1/2) In 2. But, as 
mentioned above, (|39|) is not valid below Tj> 
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The first limit corresponds to the cooling experiment with 
a sample of fixed size. By cooling it low enough we have to 
see log of some integer. The second limit corresponds to 
another set of experiments when many samples of increas- 
ing size are studied at the fixed temperature. In this case 
there is no fundamental restriction on the limiting value 
of the entropy as well as that limiting value does not have 
any fundamental meaning. 



IV. CONCLUSIONS 



Fractional entropy paradox does not exist for the finite 
system. The presented calculations provide us with the 
direct check of that. In the case of finite system Kondo 
impurity can be described in terms of free fermions. ft 
was proven above that none of them is of Majorana type. 
Using this representation we found the entropy S for ar- 
bitrary T and L. As expected, S shows a cross-over at 
Tl = v F /L. 

We derived high temperature expression for the impu- 
rity free energy and found that it coincides with the free 
energy obtained in L — oo limit. 
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